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ABSTRACT 

We present a HST/STIS spectroscopic and optical/radio imaging study of the 
Seyfert NGC 2110 aiming to measure the dynamics and understand the nature of the 
nuclear outflow in the galaxy. Previous HST studies have revealed the presence of 
a linear structure in the Narrow-Line Region (NLR) aligned with the radio jet. We 
show that this structure is strongly accelerated, probably by the jet, but is unlikely to 
be entrained in the jet flow. The ionisation properties of this structure are consistent 
with photoionisation of dusty, dense gas by the active nucleus. We present a plausible 
geometrical model for the NLR, bringing together various components of the nuclear 
environment of the galaxy. We highlight the importance of the circum-nuclear disc 
in determining the appearance of the emission line gas and the morphology of the 
jet. From the dynamics of the emission line gas, we place constraints on the acceler- 
ating mechanism of the outflow and discuss the relative importance of radio source 
synchrotron pressure, radio jet ram pressure and nuclear radiation pressure in ac- 
celerating the gas. While all three mechanisms can account for the energetics of the 
emission line gas, gravitational arguments support radio jet ram pressure as the most 
likely source of the outflow. 

Key words: galaxies: Individual: NGC 2110 — galaxies: Seyfert — galaxies: jets — 
galaxies: kinematics and dynamics — line: profiles 



1 INTRODUCTION 

Most Active galactic nuclei (AGN) are radio-quiet, i.e, 
their nuclear radio emission accounts for less than about 
5% of the bolometric luminosity of the accreting nucleus. 
However, many are still associated with non-thermal ra- 
dio sources and, in some cases, these sources can be ex- 
tended on scales of hundreds of parsecs and show col- 
limated bi-polar jet morphologies reminiscent of the jets 
seen in the more powerful class of Radio Galaxies, (e.g. 
iNagar et al.l Il999l h Among Seyfert galaxies, those that 
exhibit extended radio sources fr equently displa y accel- 
erated emission line kinematics (|Whittlel Il992l h which 
is usually interpreted as the influence of the fast mov- 
ing jet outflow on th e ionised gas in the Narrow-Line 
Region (NLR) (e.g.. 



Whittle e t all Il988l; ICapetti et all 
19961: iFalcke. Wilson fe Simpson! 1 19981: [Cooke et all l200d : 
Cecil et al.ll2002l : I Whittle fc Wilsonll2004l h 



Ionisation studies of Seyferts generally support nuclear 
UV and X-ray radiation as the principal cause of NLR 
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io nisation, even in str ongly jetted Seyferts (e.g., Mkn 78 
- I Whittle et all |2005| L Models of NLR clouds with dust 
suggest that radiation pressure from the accreting nuclear 
source plays a major part in altering the pressure equilib- 
rium of such clouds and determining their ionisation condi- 
tions dGroves et al.| [2004h In addition, the energy and mo- 
mentum deposited into the NLR from the nuclear radiation 
field can dominate the dynamics of the region and can eas- 
ily account for the energy and momentum of most outflows. 
Therefore, a key question in understanding the properties of 
AGN outflows is the relative importance of radiation force 
compared to the force from a radio jet or wind. In this paper, 
we present a detailed analysis of the kinematics and ionisa- 
tion of a nuclear outflow in the bright local Seyfert NGC 
2110 and, using geometrical and dynamical constraints, ex- 
plore and compare the accelerative capability of nuclear ra- 
diation and ram/synchrotron pressure from the jet in this 
galaxy. 

We adopt a systemic velocity of cz — 2335 km s _1 for 
NGC 2110 based on stellar absorption line measurements 



l|Nelson fc Whittle! Il995l ). With H = 72 km 
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Figur e 1. WFPC2 F606W image of NGC 2110 from iMalkarJ 
(1998). The inner 800 pc region of the galaxy studied in this 
work and displayed in Fig. 3 is indicated by a white rectangle 



this gives a distance/physical scale for the galaxy of 30.7 
Mpc/147 pc arcsec -1 . 



THE NLR OF NGC 2110: 
PREVIOUS STUDIES 



CONTEXT AND 



NGC 2110 is a nearby SO galaxy |de Vaucouleurs et al.l 
Il99ll ) with a Seyfert 2 AGN (|McClintock et al l 1 1979ft. Ini- 
tially discovered as a strong X-ray source ( Bradt et al.l 
1 19781). further spectroscopic and radio observations 
ilMcClintock et al.Ml979l ; IShuderl Il980l ; lUlvestad fc Wilson! 
Il983h revealed an extended emission line region and radio 
jet, as well as broadened lines with FWHM up to ~ 600 
km s _1 in the vicinity of the nucleus, suggesting a localised 
jet-ISM interaction. 

Compared to other Seyferts, NGC 2110 has a remark- 
ably well-defined and symmetric radio jet, with a projected 
extent of ~ 2.2" (330 pc). The inner jet is approximately 
linear on both sides of the nucleus with a north-south (PA 
0°) axis, but bends smoothly by about 40° at a projected 
distance of ~ 1", ending in large bright lobe-like structures. 
High-resolution VLBA imaging sho ws a parsec-scale featu re 
aligned with the extended inner jet |Mundell et al.ll2000l ). 

Optical continuum images show a dusty circum- 
nuclear gas disk with a well-defined spiral pattern 
l|Gonzalez Delgado et al]|2002l . also see panel 5 in Fig. 3 of 
this paper). Such disks are quite common in Seyfert galaxies 
ijRegan fc Mulchaevl[l999l ; iPogge fc Martinill2002f l and their 
spiral arms are postulated to be shocks which trigger gas 
inflow into the active nucleus. As we discuss later, the pres- 
ence of a circumnuclear gas disk may have important con- 
sequences for the morphology of the jet and the extended 
ionised gas. 

Early gro und-based studies of emission line kinematic s 
in NGC 2110 (|Wilson fc Baldwinlll985l ; IWilson et al.lll985l ) 
found evidence for rotation with the peculiar property that 
the kinematic centre was offset to the south of the contin- 
uum peak by l'.'7. This apparent paradox has been clarified 



by more recent observations. Gro und-based integral-field 
spectroscopy across the c entral 10" (jGonzalez Delgado et al.l 
120021 ; IFerruit et aT1l2004l )) has shown that the stellar veloc- 
ity field is rotationally symmetric about the optical contin- 
uum peak (the true nucleus). The gas kinematics, however, 
is confirmed to be asymm etric, though the o rigin of this 
asymmetry is still unclear. IFerruit et al.l (|2004l ) argue that 
the gas in the northern ELR is disturbed by the jet outflow, 
suppressing and lowering its average rotational velocity. 

Our study concentrates on the central 2" (300 pc) and 
uses high spatial and spectral resolution data to probe the 
properties of the inner NLR which, as we show in §4, is most 
influenced by the outflow. In §7.5, we discuss our results in 
the context of the gas flows on kpc-scales. 



3 DATASETS: REDUCTIONS AND 
MEASUREMENTS 

3.1 STIS spectroscopy 

Two high S/N STIS long-slit spectra were taken at nuclear 
and off-nuclear positions in order to cover as much of the 
central emission line region as possible. We used the G430M 
grating to the [OIIIJAA4959, 5007 doublet and H/3 at a ve- 
locity resolution of FWHM si 35 km s" 1 . See Table 1 for 
details. The [OIII]A5007 line tracks the kinematic proper- 
ties of the emission line gas, while the [OIII]/H/3 ratio is 
sensitive to its ionisation state, relatively free of reddening. 
Henceforth, we will use the labels Slit A and Slit B for the 
nuclear and off-nuclear slit respectively. 

We also use an archival STIS G 750M (red) spectrum, 
first presented in IFerruit et all [2004), which covers the lines 
of [OL1A6300, Ha, [NIL1AA6548, 6583 and [SII]AA6717, 6731. 
Like our nuclear G430M spectrum, the slit was centred on 
the continuum peak of the galaxy, though along a PA of 
155.65°. Fig. 2 compares the nuclear slit positions of the 
G430M and G750M datasets. The angular difference of 11.8° 
between the slits ensures sufficient overlap of the apertures 
across the central one arcsecond of the NLR to allow a direct 
comparison of line measurements from the red and green 
spectra in the inner NLR. 

For all the STIS data, spectrophotometric calibration 
was done using the CALSTIS pipeline and any remaining 
hot pixels were cleaned using the IRAF 'COSMICRAYS' 
task. 

Emission line properties were measured from one- 
dimensional extractions of the two-dimensional STIS spec- 
trum. The extraction widths were from one to several spatial 
increments and were designed to maximise the S/N in the 
emission lines, while isolating regions with coherent kine- 
matic behaviour. For each 1-D spectrum, we fitted a smooth 
polynomial to the underlying continuum subtracted it to 
obtain a pure emission line spectrum. The emission line 
strengths were not corrected for any underlying absorption, 
which is negligible given the very high equivalent widths of 
the lines in the spectra (> 1000 across the entire NLR). 

We measured line kinematics from the [OIII]A5007 line 
in the G430M spectra, while emission line ratios were es- 
timated using the [OIII]A5007 line as a scalable template. 
Our complete methodology is essential ly the same as that 
applied to a similar dataset for Mkn 78 l|Whittle et al.ll2005l ) 
and we refer the reader to this paper for further details. 
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Table 1. HST/STIS Spectroscopic Observations of NGC 2110 



Aperture 



Grating 

A range (A) : 

AA (A pix- 1 ) : 

Acz (km s _1 pix^ 1 ) 



G430M 

4818-5104 

0.277 

17 (@ A5007) 



NGC 2110 A Ap=52" x 0'.'2 
PA=-36.202° 
Offset=0'.'0, Of.'O 



Datasct : 
Date : 

Exposure time 



O5G401010 
24/12/2000 
1522 



NGC 2110 B Ap=52" x 0'.'2 
PA=-36.202° 
Offsct=-0'.'54, -0'.'398 



Dataset : 
Date : 

Exposure time 



O5G401020 
24/12/2000 
600 




Figure 2. The overlap between the STIS long-slit apertures 
for the G430M and G750M datasets. The apertures are plotted 
against the Ha+[NII] emission line image (see §3.2 for details). 
North to the top and East to the left. The circular inset is a zoom 
of the central 2" and certain features identified in §4 have been 
indicated. Note that the two slit apertures overlap considerably 
within a nuclear radius of ~ 0'.'5, allowing us to adequately com- 
bine line measurements from this region (the inner emission line 
plume) in our ionisation study (§5). The approximate locations 
of the six bins along each slit that were used in this study are 
shown as small white circles. 



3.2 HST and Radio Images 

Two emission line maps of NGC 2110 were prepared from 
archival HST imaging datasets. The first is a map of the 
[OIII]AA4959,5007+H/3 lines constructed from WF/PC-1 
narrow/medium-band images (Program: GO 3724, PI: Wil- 
son). The details of the instrument c onfigurations, image re - 
duction and analysis are described in lMulchaev et al.l (.1994). 
In brief, following basic calibration with the WFPC Pipeline, 
both sets of images were deconvolved with a model TINY 
TIM PSF using a Lucy- Richardson algorithm to remove the 
effects of the spherical aberration of HST. Nearby medium- 
band continuum images, also appropriately deconvolved, 
were then subtracted to get the final pure emission line im- 
ages. The plate scale of the images is 44 mas/pixel. 

A second emission line map, in Ha+[NII], was con- 
structed from archival HST images in the FR680P15 linear 



ramp filter ('on-band') and F791W filter ('off-band') taken 
on the WFPC2 PCI c amera. Calibration a nd reduction de- 
tails can be found in iFerruit et ail (|2004h . This map has 
significantly better S/N than the WF/PC-1 map and re- 
veals faint extended features in the Extended Emission Line 
Region (EELR). In addition, the F791W (approx. I-band) 
image provides an emission-line free view of the dust and 
stellar continuum geometry around the nucleus. 

For brevity, we will refer to the [OIIIJAA4959, 5007+H/3 
map simply as the '[O III]' map in the rest of this paper. 
Similarly, the H«+[NII]AA6584, 6548 map will be shortened 
to the 'Ha' map. 

A VLA X-band radio map of NGC 2110 at 8.4 GHz 
(3.6 cm), first presented in iNagar et al. I i|l999l ). was kindly 
re-reduced for our purposes by N. Nagar. The maps were 
constructed with uniform weighting and self-calibrated with 
a robustness parameter of 2, giving the best compromise 
between spatial resolution and S/N. 

3.3 Registration and Astrometry 

Accurate spatial registration of the images, radio maps and 
STIS slit positions is essential for this study. The high spatial 
resolution of HST warrants absolute positional accuracies on 
the scale of tens of milliarcseconds. This cannot be achieved 
from the HST astrometric header information alone, since 
typical absolute astrometric uncertainties in the HST focal 
plane are ~ 1". Instead, we employ a number of techniques 
to register the various datasets to each other and then an- 
chor them all to the astrometric frame of the radio maps, 
which is accurate to ~ 1 mas. 

Each STIS spectroscopic observing sequence starts with 
an acquisition image, which, for our purposes, was taken 
through the F28X50LP (long-pass) filter. The relative reg- 
istration between the slit position and the peak-up centroid 
of the acquisition image is better than 5 mas or 0.1 STIS 
pixels. The S/N of these images is sufficient for them to be 
cross-correlated with the WFPC2 F791W images, allowing 
us to match the STIS and Ha datasets together. 

We assumed that the nuclear position corresponds to 
the peak of the optical emission in the STIS acquisition im- 
ages. We independently verified this assumption by fitting 
a gaussian to the spatial profile of the stellar continuum in 
our spectra. In all cases, the peak of the continuum profile 
matched the center of the reference pixel along the cross- 
dispersion axis of the STIS spectra, as expected for a slit 
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over the peak-up centroid during the target acquisition stage 
of the STIS observations. 

To tie the STIS spectra to the WF/PC-1 images, we de- 
veloped a routine that extracted, for each STIS slit, a model 
slit from the [OIII] image and optimised its position and PA 
to match its flux profile to a gaussian-smoothed version of 
the [OIII] flux profile from the STIS spectra. The optimisa- 
tion was performed by minimising the \ 2 difference between 
the extracted and real profiles, using a robust Downhill Sim- 
plex algorithm. With a careful choice of initial parameters 
to prevent the routine from settling into local x 2 minima, 
we achieved a relative registration of the imaging and spec- 
troscopic datasets to better than 10 mas. 

In order to tie the HST datasets to the radio astro- 
metric frame, we registered the centroid of the acquisition 
image to the co-ordinates of the nuclear point source in the 
corresponding radio map, available fro m the literatur e . Thi s 
approach is justified by the result from IFerruit et all (|2004l ) 
that the kinematic center of the stellar velocity field is the 
peak in the optical light distribution. 



4 HEURISTIC DESCRIPTION 

4.1 Emission Line and Radio Morphology 

The main datasets of NGC 2110 are plotted together in 
Fig. 3. Here, we describe the principal features visible in 
the data and discuss the relationships between various iden- 
tifiable components. In what follows, we use the term 'NLR' 
to describe the entire emission line region, while the term 
EELR refers only to the outer NLR, beyond a projected nu- 
clear radius of ~ 1" (150 pc) and out to its largest extent of 
about 4" (600 pc) from the nucleus. 

The STIS spectra are sensitive enough to clearly de- 
tect the [O IIIJA5007 line within a nuclear radius of 2". 
Earlier studies have found that the kinematics of the [O 
III]-emitting gas within the inner 1" of the NLR are quite 
disturbed and turbulent, while the motion on larger scales 
is basic ally dominated by the gravitat i onal potential of th e 
galaxy l|Gonzalez Delgado et at] |2002| ; iFerruit et all l2004h . 
Therefore, probing these very central regions in detail gives 
us the best insight into the nuclear outflow. 

Within 1", both the Ha morphology (Panel 1 in Fig. 3) 
and [O III] morphology (Panel 2 in Fig. 3) share some com- 
mon features. The line-emitting gas in both images forms 
NNW a linear structure which curves over in the same di- 
rection as the radio jet and the dust lanes of the circumnu- 
clear disk (see below). The Ha image also contains a great 
deal of extended emission which is not seen in the [O III] 
image. This is primarily due to the lower S/N of the [OIII] 
image, though some of the difference is a result of the strong 
extinction towards the inner NL R south of the nucleu s from 
dust in the circumnuclear disk. IFerruit et all (120041 1 show 
that the high ionisation state of this extended gas is indica- 
tive of AGN ionisation and not distributed star formation 
in the disk. 

We turn now to the curved structure discussed above. 
Previous studies have referred to this as an emission-line 
"jet", paralleling the structure of the radio jet. The STIS 
spectra of this feature show very strong blue-shifts and com- 
plex velocity structure, which supports the notion that it is 



interacting with, or originating in, the jet. However, the in- 
ner emission line "jet" has a position angle that is offset 
by ~ 40° from the PA of the main radio jet. In fact, given 
that t he radio jet preserves its N-S direction down to parsec 
scales (|Mundell et al.ll2000h . the emission line feature prob- 
ably skirts the main radio jet, rather than being fully cospa- 
tial with it. At this juncture, we choose not to assume that 
the emission line material is accelerated by the radio jet. As 
we discuss in §8, other mechanisms may also be playing an 
important role. 

To avoid confusion with the radio jet flow, we will hence- 
forth refer to this inner emission line structure as a "plume" 
rather than a "jet". For descriptive purposes, we further 
consider the plume to have two sections: an inner and outer 
plume. The inner plume is the linear high surface bright- 
ness structure that extends 0'.'4 from the nucleus. The outer 
plume continues from the inner plume and curves to the 
NE. The motivation for separating the inner from the outer 
plume is based on their differing kinematics (§4.1). 

In order to understand the relationship between the 
radio and line-emitting components, it is important to as- 
certain the role played by dust in modifying the appear- 
ance of the emission line region. At first sight, there is a 
marked north-south asymmetry in the emission line distribu- 
tion with significantly more emission to the north. However, 
Fig. 3, dust lanes are apparent in the red continuum image 
just south of the nucleus and suggest extinction may play a 
role in creating this asymmetry. To constrain the level of ex- 
tinction, we measure the Ha/H/3 Balmer decrement south of 
the nucleus. Despite the misalignment between the G430M 
and G750M slits, there is sufficient overlap within 0'.'5 of 
the nucleus to allow a reliable estimate of Balmer decre- 
ment. Across this region, we measure Ha/H/3 — 6.1 ± 4.4, 
with large uncertainty due to the faintness of line emission in 
this region. Assuming a Case B value of 3.1, typical of AGN 
emission line regions, we estimate a V band extinction of 
Av = 1.8*i'| magnitudes, assuming the dust is distributed 
as a uniform screen. Using this Av to correct the faint [O 
III] emission south of the nucleus yields a surface bright- 
ness comparable to the northern plume, suggesting that the 
north and south inner NLR have intrinsically similar sur- 
face brightness. In reality, both the radio and emission line 
regions may have a north-south symmetry. 

Beyond the bright plume [nuclear radii > 1" (150 pc)], 
the Ha map shows curving arcs of ionised gas to the north 
and south, which are, for the most part, associated with 
dusty spiral arms visible in the continuum images. This 
EELR is significantly extended along PA ~ 160°, close to 
that of the disk major axis (and the line of nodes, assuming 
a circular disk). The k inematics of this ex tended emission 
is primarily rotational (|Ferruit et al. Il2004l ). Therefore, this 
larger-scale ionised material is probably disk gas illuminated 
by the radiation field of the central AGN, with dust and 
higher gas concentrations in the spiral arms modulating the 
appearance of the extended line emission. While the width 
of the structure is slightly narrower at the nucleus than at 
its outer extremes, its geometry does not obviously resem- 
ble a canonical do uble-sided ionisation cone found in many 
bright Seyfert 2s (|Falcke. Wilson fc Simpson] 1 1998t ). This is 
not necessarily unusual, since the illumination pattern of 
Seyfert nuclei are not always very clear. In §6, we propose a 
simple model which unifies the geometry of the EELR and 
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Figure 3. The [OIII]A5007 line for both slits, from the G430M spectra of NGC 2110. The WFPC2 Ho, WF/PC-1 [OIII] and WFPC2 
F606W images are shown for reference. The slit positions are indicated. Radio contours are over-plotted on the [OIII] image. The 
hatched ellipse represents the FHWM beam of the radio map. Dotted lines (dashed for the nucleus) are drawn at 1" intervals along the 
slit direction. 
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the radio jet, and reconciles the illumination pattern with a 
bi-cone. 



4.2 Emission Line Kinematics 

In Fig. 4, we plot peak velocity with respect to systemic 
(Panel 1) and FWHM (Panel 2) measured from the [O 
III]A5007 line, as a function of nuclear distance along both 
slits. When comparing these measurements to the panels of 
Fig. 3, clear relationships can be seen between the kinemat- 
ics of the line-emitting gas and its spatial structure. A com- 
pact knot of gas surrounds the nucleus in both [O III] and 
Ha maps, resolved with HST, but only marginally resolved 
by the STIS spectrograph due to the Of.' 2 slitwidth. This 
emission line knot displays extremely broad lines (FWHM 
of 1100 km s _1 ), a median biueshift of 250 km s _1 and 
a predominantly blue wing. In addition, in spatial incre- 
ments south of the nucleus, a substantial red wing to the 
[O III] profile becomes evident. This is consistent with an 
unresolved bipolar nuclear outflow with a dust-extincted red 
component, smeared out by the instrumental PSF. 

What might be the origin of the very broad central 
component? iMoran et alj (|2007l ) estimate a central super- 
massive black hole mass in NGC 2110 of M B h = 2 x 10 s 
Mq from its central velocity dispersion of <r» = 220 km 
s -1 jNelson Sz Whittle! [l995). using the relationship from 
iTremaine et al] ] 20021 ). The radius of gravitational influence 
for a black hole of this mass given by: 



GMbh 



(1) 



which is 18 pc, or 2.4 STIS spatial increments. While ap- 
proximate and dependent on the shape of the potential, this 
estimate of R g closely matches the size of the nuclear com- 
pact knot. It is quite possible that some of the high velocities 



measured in this knot can be attributed to gas within a few 
parsecs of the black hole. On the other hand, there appears 
to be a continuity in the bulk velocity of the gas across the 
transition between the compact knot and the rest of the 
inner plume, implying that the plume is the resolved exten- 
sion of the same bipolar outflow that extends down into the 
knot. Also, as discussed in §5.3, the ionisation properties of 
the knot are similar to that of the inner plume. Therefore, 
the knot probably contains gas in rotation around the black 
hole as well as an outflowing component. We discuss the 
outflow and its accelerating mechanism in §7. 

Beyond the core, the blueshifted gas extends along the 
entire length of the inner plume, decelerating with nuclear 
distance. The FWHM of the gas drops from about 900 km 
s _1 within O'.'l to systemic at the end of the plume where the 
bulk blue shift drops to zero. Past the end of the inner plume, 
the linewidths rapidly become narrow and the gas returns to 
a quiescent state, consistent with almost normal rotation in 
the circumnuclear disk. This velocity pattern suggests that 
the inner plume is a site of strongly outflowing gas. Note that 
the inner plume lies well outside the sphere of the influence 
of the black hole, so its rapid deceleration should not be 
interpreted as a Keplerian drop-off. 

What about the outer plume? Does it inherit the dis- 
turbed motion of the inner plume? Slit B samples the outer 
plume as it intersects the central axis of the radio jet. This 
is, in fact, the only part of the plume structure that appears, 
in projection, entirely within the confines of the radio jet. 
Surprisingly, the ionised material here has narrow linewidths 
and appears to be quiescent and rotational, similar to the 
gas beyond the plume in Slit A. In contrast to the inner 
plume, the outer plume is likely to be an AGN-illuminated 
nuclear spiral arm, lying either in front or behind the jet. 
The continuum image (Panel 5 in Fig. 3) supports this idea: 
unlike the inner plume, the outer plume is associated with 
a similarly shaped dust feature. As we discuss in §6, an in- 
vestigation into the relative geometry of the NLR and the 
radio jet provides further evidence that the jet is viewed in 
projection against the outer plume. 

Our preliminary analysis of the interaction indicates 
that the inner plume is highly disturbed and directly ac- 
celerated by the nuclear outflow. In the next section, we 
consider the physical nature and excitation of the ionised 
gas in the nuclear knot and the inner plume. 



5 IONISATION CONDITIONS 
5.1 Ionisation models 

A knowledge of the ionisation conditions of the emission 
line gas helps to set the stage for a study of the energet- 
ics of the inner NLR. The radiation field of the central 
AGN is often a major source of thermal and kinetic energy, 
through photoionisation heating and radiation pressure. In 
addition, strong shocks, driven into the NLR clouds by inter- 
actions with the radio jet o r outflow, can also be a significant 
source of ionising photons dViegas-A ldrovan di fc Gruenwaldl 
ll98Sl ; lDopita fc Sutherlandlll996l : lAllen et al.ll200Sl ) and. urT 
der certain conditions, this can dominate over nuclear pho- 
toionisation. Since shocks are sinks of the jet kinetic energy 
flux, constraints from shock ionisation can be used to deter- 
mine jet energies. In addition, variations in ionisation state 
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Figure 5. Emission line ratio-ratio plots for the inner plume in NGC 2110. Data points are plotted as filled circles with 2a error bars. 
Two sets of ionisation model tracks are plotted: U,; ust models (dashed lines - sequences in ionisation parameter) and shock models from 
lAllen et all ^0081 1 (solid lines - sequences in shock velocity). Symbols plotted on each sequence help identify the parameter value. Shock 
velocities range from VJ,j, = 400 km s" 1 (cross) to 1000 km s _1 (square), with plus points every 100 km s — 1 and a star at 700 km s _1 . 
\ogUdust ranges from —4.0 (cross) to 0.0 (warped square), with a plus point every 0.5 dex and a star at —2.5. Thin/thick lines are pure 
shock/shock+prccursor models for the shock sequences and 2x/4x solar for the Udust sequences. More details can be found in §5.1. 



can arise from changes in gas density and Lyman continuum 
optical depth. A comparison of line strengths with the pre- 
dictions of models can throw light on the micro-scale nature 
of NLR gas, which then feeds into a dynamical analysis and 
estimates of jet properties. 

We employ two grids of ionisation models applicable 
to AGN emission line regions. The first considers the pho- 
toionisation of dusty, optically-thick gas by the EUV and 



X-ray contin uum of the n u clear source. These models were 
presented in iGroves et al.l (|2004 ) and span a range in NLR 
metallicity, mean ionisation parameter U, gas density and 
ionising continuum slope. We call these the Udust models. 
The second set of models are those of autoionising radia- 
tive shocks, i.e, shocks that ionise an optically thin pre- 
cu rsor region o f unsh ocked g ELS, ELS presented and discussed 
in Iau en et alj (^OOct ). These models extend the work of 
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iDopita fc Sutherland! (|l996h . spanning a greater range of 
shock velocity, gas density, gas metallicity and magnetic 
compression parameter. Our approach is to compare mea- 
sured line ratios with the predictions of both of these models, 
in terms of their location on line ratio-ratio diagrams and, 
since we have spatially resolved the inner NLR, trends in the 
measurements. A comparison of several different diagrams 
can allow an assessment of the importance of these two very 
different ionisation mechanisms. 



5.2 Line ratio-ratio diagrams 

Fig. 5 plots combinations of optical emission line ratios mea- 
sured for the inner plume. It is worth recalling that these di- 
agrams compare line ratios from two different STIS datasets, 
though any given ratio is constructed from lines measured 
in a single dataset (G430M or G750M). The misaligned slits 
(angular offset of about 12°) sample slightly different parts 
of the inner plume, from almost full overlap at the nucleus, to 
approximately 25% overlap at the edge of the plume (Fig. 2). 
Our approach will be valid as long as there are no strong ion- 
isation gradient perpendicular to the plume, which seems 
very unlikely. 

Line ratios, calculated from the flux measurements in 
Table 2, are plotted in Fig. 5 as filled circles with 2a error 
bars. In what follows, we exclude the set of STIS spatial in- 
crements on or immediately adjacent to the nucleus to avoid 
complications from the smearing of the nuclear point source 
due to the 0'.'2 width of the slit. The approximate central lo- 
cations of these increments are plotted in Fig. 2. As a guide 
to the reader, the data points with lowest [OIII] A5007/H/3 
are the closest to the nucleus. 

We have compared the full space of parameters from 
both Udust and shock models against many combinations 
of line ratios. Fig. 5 is a distillation of a set of 4 line ra- 
tio diagrams and a small subset of model tracks that bring 
out our main results. Plotted using dashed lines, the rep- 
resentative Udust models feature a power-law ionising spec- 
trum in frequency, with a spectral index of —1.2 (typical of 
Seyferts), illuminating gas clouds with an electron density 
of 1000 cm -3 , with tracks in ionisation parameter spanning 
—4.0 < log U < 0.0 and metallicities of twice and four times 
solar. Representative pure shock and shock+precursor mod- 
els are plotted as solid lines, with pre-shock densities of 1 
cm -3 and shock velocities between 400 and 1000 km s -1 
as well as high magnetic fields 10 mG (highly magnetised 
shocks). We find that in order to match the strengths of the 
[S II] and [N II] lines, we require super solar abundances in 
both model sets. These models were chosen to span the mea- 
sured line ratios, while fixing parameters that have a minor 
influence on the position of the model tracks on these dia- 
grams. We refer the reader to the original papers for a full 
description of the parameter space covered by the models. 

A examination of the line ratio-ratio diagrams yields 
the following conclusions: 

a.) A gradient in the ionisation sensitive line ratio [O 
III]/H/3 is observed, towards lower ionisation at smaller 
radii. The low ionisation ratios of [N II] and [S II] with Ha 
show no trend and indeed, the entire plume occupies a very 
small region of the diagram that involves combinations of 
these three lines (Fig. 5b). The [O I] /Ha line ratio exhibits 



an anti-correlation with [O III]/H/3, which, interestingly, is 
not mirrored by any of the model tracks. 

b.) Considering the shock models, the low ionisation 
points prefer pure shocks with shock velocities around 
700 km s _1 while the high ionisation points match 
shock+precursor models with higher velocities. The tracks 
for unmagnetised shocks do not span the range of line ra- 
tios as well - weakly magnetic pure shocks produce too weak 
[O III]/H/3, while weakly magnetic shock+precursor models 
produce low values of [S II], [N II] and [O I] relative to Ha. 
Note that the transition from having a significant precur- 
sor contribution near the nucleus to pure shocks at larger 
radii is required to match the trends in the measurements. 
Neither a full set of pure shock or shock+precursor mod- 
els actually spans the gap between low and high ionisation 
points. Mod els with hig h er ab undances, while not explicitly 
modeled bv lAllen et al.l |2009), are unlikely to improve the 
match to the measurements, as they will move the [O III] 
strengths higher in the shock+precursor models and away 
from the measurements. 

c) On the other hand, the dusty photoionisation mod- 
els can consistently reproduce the measured ratios with 
\ogU ~ 2.0 and abundances about 3 x solar. There may be 
a preference for a harder ionising power-law than we have 
adopted, which will boost [S II] and [N II] relative to Ha, 
yielding a better fit in Fig. 5b. A harder ionising contin- 
uum will also lead to a better match to t he [O I] strengths 
in Fig. 5d. However, iGroves et al.l (|2004l 'l do not consider 
power-law indices smaller than —1.2 and modelling this sce- 
nario is beyond the scope of this paper. In the context of 
the photoionisation, the gradient in [O III]/H/3 implies lower 
ionisation parameters at smaller nuclear distances. Since the 
ionising field is expected to decrease with distance, due to 
geometric dilution and increasing net optical depth, the gra- 
dient implies that the density of the ionised gas increases 
sharply towards the nucleus. 

For a galaxy with the luminosity of NGC 2110 
(Mb = —19 . 3), lo cal metallicity scaling relationships (e.g., 
ISalzer et al.l I2005T 1 predict oxygen abundances around so- 
lar. The inner kpc of NGC 2110 appears to have signifi- 
cantly higher metalllicities than the mean galaxy popula- 
tion, though an excess of 0.5 dex is not much larger than la 
scatter in abundance for galaxies of this luminosity. 

d) Unlike the other ratios, [O I] /Ha is not fully con- 
sistent with the best parameters for the Udust sequences 
(Fig. 5d). Shocks give a somewhat better match, though, 
as before, the gradient supports a switch from pure shocks 
to shock+precursor emission over the length of the plume. 
However, [O I] is a very difficult line to model because its rel- 
ative strength is influenced considerably by the high energy 
(soft X-ray) contribution to the ionising spectrum, which 
produce the partially ionised zones needed to excite neu- 
tral oxygen significantly. NGC 2110 has an extended soft X- 
ray excess compar ed to normal Seyferts (|Weaver et al.|[l995l ; 
lEvans et al.ll2006h and this may explain why [O I] is stronger 
than predicted by the photoionisation models. 

5.3 Global Constraints 

A simple check on the feasibility of nuclear photoionisation 
can be made by comparing the mid/far IR luminosity of 
the galaxy with the ionising luminosity necessary to power 
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Table 2. Emission Line Flux Measurements of the the Inner Plume 



Nuc. Dist. All fluxes are in units of 10 13 erg s 1 cm 2 arcsec 2 



(arcsec) H/3 [O III]A5Q07 [O I]A6300 Ha [N IIJA6584 [S II]A6717 + 6731 



0.08 


0.79 ±0.24 


3.82 ±0.23 


2.27 ±0.05 


3.50 ±0.52 


5.71 ±0.57 


4.02 ± 1.01 


0.13 


0.52 ±0.21 


3.73 ±0.22 


1.14 ±0.05 


2.13 ± 0.32 


3.50 ±0.35 


2.66 ± 0.40 


0.18 


0.40 ±0.16 


3.76 ±0.23 


0.79 ± 0.03 


1.66 ± 0.25 


2.67 ±0.40 


1.85 ±0.28 


0.23 


0.34 ±0.14 


3.61 ±0.18 


0.67 ±0.03 


1.61 ± 0.24 


2.50 ±0.25 


1.58 ±0.32 


0.28 


0.38 ±0.15 


3.87 ±0.19 


0.61 ±0.04 


1.41 ± 0.21 


2.42 ± 0.24 


1.46 ±0.29 


0.33 


0.38 ±0.13 


4.10 ±0.12 


0.67 ±0.05 


1.43 ± 0.21 


2.61 ± 0.26 


1.53 ±0.28 



the observed line emission. Seyferts are known to be power- 
ful IR sources with particularly warm col ours l|Riekelll97g| ; 
iPerez-Garcfa fc Rodrfguez-Espinosal 120011 ). Most of the IR 
emission is believed to come from the AGN's optical-UV- 
X-ray output that is intercepted and reprocessed by warm 
dust near the nucleus and in the torus. Since NGC 2110 does 
not show any signs of strong star formation, either circum- 
nuclear or in the large-scale galaxy disk, its IR luminosity 
Lib. (8 — 1000/im) should be do minated by the output of t he 
AGN. Using the calibrations of lSanders fc Mirabej (|l996h : 



Lir « 1.8 x KT 11 4vrd 2 
x (13.5Sia + 5.2S 2S + 2.585 , 6o + S100) erg s" 1 (2) 



where S12 = 0.35, S25 = 0.84, S 60 = 4.13 and S100 = 5.68, 
all in Jy, are the fluxes of NGC 2110 in the four IRAS 
bands and d cm is the distance to the galaxy. This gives us 
log Lm ~ 43.7 erg s~ . Assuming that most of the AGN's 
radiation comes out at UV to X-ray wavelengths and is re- 
processed, its bolometric luminosity L p h ~ Lie,- This esti- 
mate is consistent with results from spectral fitting of the 
nuclear X-ray emi ssion, appropriately scaled using a b olo- 
metric correction (| Evans et al. l l2006l ; lMoran et al.ll2007l ). 

One can estimate the total ionising luminosity using the 
total H/3 flux of the Seyfert. Assuming a power-law AGN 
spectrum of the form L„ <x v~ a ° (a > 1.0), a dust absorp- 
tion cut-off at Xi ow and a covering factor Q, d for the ionised 



JH/3 



^1 Ci L 



ph ~ 



OtH/3 t H p a 



OiB tH 



o-l ( 912 Y"- 1 



0.023 n d L ph — (t— ergs" 1 (3) 



where olhp and as are the H/3 and total case B recombina- 
tion coefficients for hydrogen at T — 10 4 K, and eup = hvup 
is the average energy of an H/3 photon. 

Q d is the fraction of 4ir steradians intercepted by line 
emitting gas. In ^3 we calculate the intrinsic cloud covering 
factor Ci ~ 1. Taking the ionising cone half-angle to be 45° 
(based on estimates in §6), we calculate ild ~ 0.3. With 
canonical values of a = 1.2, Xi ow = 4000A and using our 
estimates of L p h and QCi, log Lap ~ 40.8 erg s _1 . Despite 
the uncertainties in this approach, this estimate compares 
very well with the total extincti on-corrected H/3 lum inosity 
of logLiip ~ 40.8 erg s" 1 from IWilson et all l| 19851 ) based 
on early ground-based long-slit work. 

The relative importance of shock ionisation can be 



tested in a similar fashion. The Extended Emission-Line Re- 
gion (EELR: nuclear radii > 2") is undoubtedly ionised by 
the central AGN, since the undisturbed gas kinematics on 
these scales cannot sustain the fast shock velocities (> 500 
km s" 1 ) needed to produce its high level of excitation. 

We concentrate instead on the disturbed inner plume. 
We can compare the ionising flux experienced by gas in 
the plume both from shocks and the central source. We 
can estimate the amount of ionising UV radiation gener- 
ated per unit shock area, u sing the following relation from 
iDopita fc Sutherland! (|l996l ): 



fuv,shock = 1.11 x 10 3 n Vsh 4, erg cm 2 s 1 



(4) 



where V s h is the shock velocity in units of 100 km s -1 and 
n is the preshock density in cm -3 . For V s h = 500 km s _1 
and n = 1, fuv,shock = 0.14 erg cm" 2 s -1 . 

Taking the photon luminosity of the AGN L p h and a 
characteristic distance of the inner plume (0'.'5 « 75 pc), we 
can estimate an ionising photon flux per unit area from the 
nuclear source at the plume of fuv,nuc ~ 75 erg cm -2 s _1 , 
almost 3 orders of magnitude higher than the UV flux from 
shocks. This suggest that the ionisation of the inner plume 
is dominated by the photon output of the nucleus. Taken 
together with the dominance of nuclear photoionisation in 
the EELR, we conclude that shock ionisation contributes 
very little to the line emission in NGC 2110. 



5.4 Line Profiles 

A different handle on the nature of the ionisation comes from 
a comparison of line profiles in the emission line plume. In 
Fig. 6 we plot the normalised profiles of [O III]A5007 and H/3 
integrated over the inner plume, the nuclear knot and the 
region of weak emission south-east of the nucleus along Slit 
A. A measure of the ionisation level as a function of velocity 
is revealed by the relative strengths of [O III] to H/J across 
the profiles. 

In all three regions, the [O III] line is wider than H/3. In 
the nuclear knot and inner plume, this difference in linewidth 
takes the form of a stronger, more highly ionised blue wing 
on [O III]. This blue wing in the nuclear knot forms the 
base of the decelerating inner plume emission and displays 
a similar line ratio ([O III]/H/3~ 4.0), implying a continuity 
between these regions. The line core, with absolute velocities 
< 200 km s" 1 around systemic, exhibits a low LINER-like 
ratio ([O III]/H/3~ 2.8). In contrast to the inner plume, the 
red half of the [O III] line in the south-east region may be 
stronger than H/3 though the low S/N in this region pre- 
cludes a definite measurement. This suggests that a highly 
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extincted plume-like structure also exists in the inner south- 
east NLR, consistent with the results of §4.1. 

Based on these trends, the following picture presents 
itself: the low velocity line core is emitted by ionised gas 
that is relatively undisturbed and might lie in the plane of 
the circumnuclear disk. The outflow, in contrast, is highly 
ionised and possibly mass loaded from material accelerated 
off the low ionisation clouds in the disk. This enhancement 
of ionisation could be due to a greater fraction of low den- 
sity, optically thin gas in the outflow. As it expands away 
from the nucleus, the outflow decelerates rapidly and finally 
merges with the rest of the disk, losing its distinct identity 
and kinematics. 

To conclude, the line-ratio analysis and the agreement 
between the line and reprocessed ionising luminosities of the 
AGN lends some support for widespread nuclear photoioni- 
sation. In addition, the presence of strong trends in ionisa- 
tion across line profiles in the innermost plume allow us to 
differentiate between an undisturbed, LINER-like disk com- 
ponent photoionised by the AGN, and a more highly ionised 
outflow component. 



6 DESCRIPTION OF THE INTERACTION 

A geometric and kinematic model for the NLR is necessary 
for a dynamical analysis of the interaction. To this end, we 
first introduce various structural components of the NLR 
and then discuss how they relate to each other. 

Based o n elliptical isophote fits to the F791W contin- 
uum image, iGonzalez Delgado et all l|2002l ) estimate a cir- 
cumnuclear disk inclination of approximately 42°, with the 
line of nodes aligned along PA 163? 5. The contrast of the 
dust lanes indicates that the western side of the disk is in- 
clined towards the line of sight. Combining this information 
with the direction of rotation of the ionised gas in our spec- 
tra (blueshifted to the north and redshifted to the south), 
we deduce that the spiral arms trail the disk rotation in the 
normal fashion of large scale galaxy spirals. 

From Fig. 3, a comparison of the Ha map (Panel 1) and 
the radio map (contours in Panel 2) shows that the major 
axis of the ionised gas distribution is inclined to the axis 
of the radio jet by about 20°. At face value, this implies 
a misalignment between the radio jet ejection axis and the 
normal to the plane of the putative central molecular torus, 
which defines the orientation of the ionisation cones. While 
this can be explained by models with accretion disk pre- 
cession or warps, we propose the following geometric model 
to account for the misalignment (see Fig. 7). We begin by 
assuming that the inner jet is in fact aligned with the axis 
of an illumination cone. The cone axis is inclined from the 
normal to the circumnuclear gas disk by an angle <j>. In the 
assumption that almost all the circumnuclear gas in con- 
fined to the disk and that the disk scale-height is small (a 
thin disk), the emission line distribution is determined by 
the ionised gas in the disk, which takes the form of a planar 
section through the illumination cone (a conic, though only 
approximately, given the inhomogeneity and patchy dust ob- 
scuration within the disk) . The axis of the EELR is then the 
projection of the illumination cone axis (i.e., the inner radio 
jet axis) onto the circumnuclear disk, with the direction of 
projection along the normal to the disk plane. 
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Figure 6. A comparison of the [O IIIJA5007 and H/3 line profiles 
in three regions of the inner NLR: the inner plume (top), the 
nuclear knot (centre) and the dust obscured region just south-east 
of the nuclear knot along Slit A (bottom). Note the distinctive 
trend in the upper two panels for the blue wings of the profile to 
be stronger in [O III] compared to H/3, implying a high degree of 
ionisation. A hint of the reverse (a higher ionisation red wing) is 
visible in the south-east inner NLR, but the S /N of the lines is 
very low here, due to strong dust extinction. 

We define two mutually rotated reference frames: a) a 
'disk frame' with the circumnuclear disk in the XY plane 
and the normal to the disk along the Z-axis, and b) a 'sky 
frame', with the sky in the X'Y' plane and the line of sight 
to the nucleus along the Z'-axis. We choose the X and X' 
axes to lie along the line of nodes of the disk, making it 
the common axis of rotation between the two frames. We 
denote the rotation angle (anti-clockwise looking towards 
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the origin) to be a. This angle is the inclination angle of the 
disk . The jet axis is taken to be a radial vector with an angle 
from the Z-axis of <j> (<j> ) and from the line of nodes of 9 (9 ) 
in the disk (sky) frames. Basic rotational transformations 
between spherical co-ordinate frames gives us relationships 
between these angles: 



tan f 



cos a sin 9 sin < 



sm a cos < 



cos H sm <p 
cos a cos d> + sin a sin 9 sin < 



(5) 
(6) 



Since the axis of the EELR (i.e., the projection of the 
jet axis onto the XY plane in the disk frame) is essentially 
aligned with the line of nodes, then 9 « 0° (see §4.1). Taking 
a — 42° (the inclination of the circumnuclear disk) and 9 ~ 
— 16? 5 (the angle between the radio jet axis and the line of 
nodes in the sky frame, which is equal to the PA of the line 
of nodes, since the jet is aligned north-south), we can solve 
for (f> using Eqns. 5 and 6: 



be excluded, it seems unlik ely given the constan t jet align- 
ment down to parsec scales (|Mundell et alJl200Ch . Given the 
strong gradient in gas density perpendicular to the circum- 
nuclear disk and the fact that the jet bends towards the 
minor axis of the disk, we propose that the jet bends as a 
result of refraction in the vertical pressure gradient of the 
disk. 

With this configuration for the jet and the disk, the 
inner emission line plume marks the location of higher den- 
sity outflowing gas confined to the disk, ionised mostly by 
UV radiation from the AGN. We find little evidence for gas 
entrained into the jet itself: for most of its extent, the jet 
and its terminal lobes are free of associated high velocity 
emission line material. This may be compared to the view 
of the jet interaction in Mkn 78 (| Whittle fc Wilson! [20041 '). 
where large amounts of line-emitting gas are embedded and 
carried along by the jet flow. With this simple model, we 
now explore the general energetics of the NLR and radio jet 
from a dynamical standpoint. 



tan < 



cos a sin 9 — cos 9 tan 9' tan ( 



which gives <f> as 70°. In other words, the jet axis is inclined 
to the plane of the disk by 90° - <j> « 20° Using Eqn. 6, 
we calculate the angle between the jet axis and line-of-sight 
4 ~ 72°. 

In the limit of a very thin disk and a well-defined illumi- 
nation cone, the opening angle of the cone can be estimated 
from the apparent opening angle of the ionised section of the 
disk. Taking this to be between 70° to 90° , simple geometric 
arguments give cone opening angles of around 80° to 100°, 
typical of Seyfert ionisation cones. 

Finally, an estimate can be made of the thickness of the 
disk. If we assume that the kinematically disturbed inner 
plume coincides with the projected length I of the outflow 
as it interacts with the disk, then the disk thickness D is 
given by: 



D = 



21 



150 pc 



(8) 



with / = 0'.'5 (75 pc) and the values of <j> and <j> estimated 
above. This number is quite plausible, though some caveats 
must be noted. Real disks are not sharply edged but likely 
to have an exponential vertical gradient in gas density, while 
the gas distribution in the disk will be highly clumpy and 
non-uniform. Both these effects would lead us to underesti- 
mate the disk thickness. Another important consideration is 
that the outflow could disrupt the orderly structure of the 
cold disk gas. 

The nature of the circumnuclear disk has also im- 
plications for the radi o source morphology as well. 
lUlvestad fc Wilson! l| 19831 ) consider possible causes for the 
symmetric bending of the jet in NGC 2110: a) ram pressure 
forces from the rotating ISM, b) static pressure gradients 
in the bulge and, c) precession of the jet axis. It is possible 
to rule out ISM ram pressure as the main cause of the jet 
bend, since the inferred direction of the disk rotation is op- 
posite to that required to bend the jet to its current shape, 
except for the unlikely case where the tangent plane to the 
jet bea m sweeps across the LOS plane over the extent of 
the jet (jFiedler fc Henriksenlll984h . While precession cannot 



7 DYNAMICAL STUDY 

NGC 2110 hosts a significant radio jet and a relatively lumi- 
nous AGN, both of which can drive a fast moving outflow, 
either through the action of relativistic pressure, ram pres- 
sure or the momentum imparted by nuclear radiation, i.e, 
radiation pressure. In this section, we estimate the capacity 
of each of these mechanisms to provide the total energy and 
momentum of the ionised outflow, in order to gauge their 
importance. 



7.1 Emission Line Region Properties 

If the gas is primarily photoionised (§5), the masses, mo- 
menta and energies of the inner plu me can be derived from 
basic nebular photoionisation theory (|Osterbrockl 1989T ) . The 
total flux in H/3 measured from the STIS spectra and cor- 
rected for an extinction of Av = 1.1 magnitudes is Fap ~ 
2.4 x 10 -14 erg cm" 2 s -1 . From this, the mass in emission 
line gas can be estimated using the following relation: 

M em » 171 F H p d 2 (9) 

giving M em ~ 1.5 x 10 4 solar masses, where we have 
adopted a hydrogen number density of n em = 1200 cm -3 , 
derived from the r atio of the [S II] AA6716, 6731 doublet 
jFerruit et al.|[l999l ). Approximating the plume as a cuboid 
of volume V em « 0.2 x 0.2 x 0.5 arcsec 3 (or 6.76 x 10 4 pc 3 ), 
we can estimate the filling factor of ionised clouds: 



ffe 



TFlp Vem 



(10) 



to be 8 x 10 -3 , a typical value for NLR clouds. Combining 
//em with an estimate for the column density of the NLR 
clouds Nem ~ 2 x 1 Q 21 cm -2 from X-ray absorption models 



|Ferruit et all fl999). we derive typical covering factors of 
ionised gas in the NLR: 



3 TLem I 



2 N e . 



(11) 



In addition, a simple estimate can be made of the age 
of the interaction from the approximate time taken by the 
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Figure 7. A schematic diagram of the circumnuclear region of NGC2110. A inclined rotating disk containing gas and dust is intersected 
by a biconical radiation field. A radio jet emerges along the axis of the bi-cone but bends away from the disk. The various angles and 
axes are defined in §6. 



emission line gas to cross the region of the plume. With a 
characteristic velocity v em « 400 km s _1 , we get a timescale 
t ag e ~ l/v em ~ 2 x 10 5 years. This age estimate is essentially 
an upper limit if the gas in the plume is decelerating with 
distance from the nucleus - a plausible interpretation of the 
bulk velocity profile of the plume. 

Combining M em with the kinematics of the [O III] line 
and accounting for a projection angle to the line of sight 
of (j) = 70°, the total (translational and internal) emission 
line kinetic energy is E ke ~ 10 52 - 3 erg and a bulk mo- 
mentum is G em ~ 10 44 ' 8 dyne-s. Ek e may be compared 
to the equipartition energy in the radio source (~ 10 53 



erg. lUlvestad fc Wil son 1983) or the integrated energy from 
nuclear UV radiation over the lifetime of the interaction 
(~ L p h x tage = 5 x 10 56 erg). Energetically, both the jet 
and nuclear radiation are capable of supplying the energy of 
the outflow, with the nuclear radiation field dominating, by 
far, the energy budget of the system. 

Assuming that the relativistic particles and magnetic 
fields of the radio-emitting plasma are in equipartition, var- 
ious physical param eters o f the r adio jet can be estimated 
using relations from iMileyl dl98Ch and radio measurements 
from lUlvestad fc Wilson! (|l983T ). Radio source pressures are 
modest, with log P re i ~ —9.0 dyne cm -2 . These may be 
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compared to the pressures of the emission line gas and the 
extended hot gas in the NLR responsible for the soft X-ray 
emission. Using n em and a temp erature T em ~ 1360 K from 
the [O III] AA4363/5007 ratio (|Ferruit et alj|l999h . we de- 
termi ne the ionised gas p ressure to be log P em ~ —8.6 dyne 
cm" 2 . [Evans et all (|2006T ) estimate the X-ray emitting gas 
to be at a pressure of log Px ~ —8.5 Both the emission line 
gas and hot gas appear to be slightly over-pressured with 
respect to the radio jet plasma, but, given the uncertainties 
of a few associated with such calc ulations, this is basically in 
agreement with the conclusion of lUlvestad fc Wilson! (|l983T l 
that the phases are in rough pressure balance. Note that 
there is no good evidence that equipartition is a valid con- 
dition in such weak Seyfert jets. Studies of FRI jets in radio 
galaxies suggest that they are particle-dominated, probably 
due to entrainment of a heavy t hermal ions, and quite out o f 
equipartition (see the review bv lMcNamara fc Nulsenll200it ). 
It is unclear whether this scenario will also apply to small 
Seyfert jets in gas-rich circumnuclear environments. How- 
ever, given the similarity of the pressures across phases in 
NGC 2110, we assume that equipartition is probably valid 
and proceed accordingly. 

7.2 Possible Acceleration Mechanisms 

Over the lifetime of the outflow, the forces that accelerate 
the outflow have to be able to supply the bulk momentum 
measured in the inner plume. The momentum supplied to 
the gas per unit time, the momentum flux H em = G em /t ag e, 
must be matched by any relevant accelerating mechanism. 
From the estimates above, logil em w 32 dyne. 

A. ) Relativistic Pressure 

Is the radio source pressure capable of providing the mo- 
mentum of the emission line gas? We envision the jet pres- 
sure acting on a sheath of ionised gas in contact with the jet 
surface and accelerating the gas outward from the surface. 
The area of this sheath Abeam is approximately the surface 
area of a cylindrical beam of length I and cross-sectional ra- 
dius of the jet Rj « 70 pc, measured from the radio map. 
The relativistic pressure acting over this area gives a mo- 
mentum flux Of n re j m PrelAbeam, OT log Yl r el ~ 31.8 dyne. 

Within the considerable uncertainties, H re i ~ n em and the 
relativistic pressure could, in principle, accelerate the emis- 
sion line gas. 

B. ) Jet Ram Pressure 

The ram pressure of the fast moving jet flow as it prop- 
agates from the nucleus will accelerate clouds that lie along 
its path. While the physics of interactions between a jet and 
ISM clouds involves complex hydrodynamic processes such 
as shocks, cloud-breakup and ablation, we take a broader 
view and ask simply whether a jet with a ram pressure 
Pram = PjVf can satisfy the constraints of IT em , where pj 
is the jet density and Vj is the jet velocity. 

A constraint on these jet parameters comes from the 
morphology of the jet. We have argued that the most likely 
explanation for the symmetrical bend in the jet is the action 
of a disk pressure gradient. Modelling the disk as a slab 
with a scale height much smaller than the disk plane scale 
length, we can appl y the following condit ion for significant 
jet deflection from lHenriksen et al.l (| 198lh : 



P> V ? * ^J Ps (12) 

where P s is the disk pressure at the sonic point of the jet and 
7 is the polytropic equation of state of the jet material (= 
4/3 for a relativistic gas). Assuming that the line emitting 
gas is in rough pressure equilibrium with the hot atmosphere 
of the galaxy, i.e, P em ~ P s , we get jet ram pressures of log 
Pram ~ —8.0 dyne cm -2 . This pressure acting over the area 
of cross-section of the jet gives a force of II ram « 7r7? 2 P ram 
or log Uram ~ 33 dyne. Interestingly, requiring the jet to 
bend in the disk pressure gradient also implies a jet velocity 
and density that can easily accelerate the outflow by ram 
pressure. 

Taking a fiducial proton number density for the jet of 
Pj/rrip — 10 -2 cm -3 , we get Vj ~ 10 _2 c, which is highly sub- 
relativistic. In other words, the requirement of substantial 
buoyant force on the jet material in the galaxy's atmosphere 
implies that the jet is either highly underdense or highly 
subrelativistic. 

C.) Radiation Pressure 

The emission line gas in the plume is also illuminated 
by a powerful ionising radiation field from the nucleus. The 
pressure from this radiation is potentially a major source 
of momentum to the ionised gas. We estimate its momen- 
tum flux from the photon luminosity of the nucleus and the 
covering factor of the inner plume fid ~ 0.02: 

I7w = (13) 

c 

which yields logII ra[ j ~ 31.5, a little lower than H re i- Within 
our uncertainties, radiation pressure is also capable of driv- 
ing the nuclear outflow. However, the ability of radiation 
pressure to impart momentum to the ionised gas depends 
on nuclear distance simply because at smaller distances the 
covering factor fid is greater. 



7.3 The Mass Outflow Rate 

The mass outflow rate in the jet is given by: 

Mjet = K Pj VjRj = nPramlfj /Vj w 7.8xlO _4 /Vj M Q yr 

where we have adopted the value of Pram from the jet bend 
analysis above. 

The Eddington accretion rate for the SMBH in NGC 
2110 is ~ 5.3 MQyr" 1 , assuming a radiative efficiency of 0.1, 
while the mass accretion rate for the estimated bolometric 
luminosity of the AGN (L ph ) is M acc ~ 8.8 x 10" 3 M yr _1 . 
Interestingly, the mass accretion rate of the SMBH is close to 
(possibly even lower than) the mass outflow rate in the jet if 
we adopt our fiducial value for Vj . As long as the constraint 
on Pram is valid, this may suggest that the jet carries away 
a substantial fraction of the mass accreted onto the SMBH. 
However, it is also possible that the current AGN luminosity 
may have varied substantially over tag?, since the time of 
the ejection of the jet, that the jet is highly underdense or 
modestly relativistic, or that the jet ram pressure estimate is 
too high. The large number of contributing factors makes the 
determination of the most likely cause for this discrepancy 
difficult to ascertain. 
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7.4 The Role of Gravity 

If any of the acceleration mechanisms discussed above are 
to accelerate gas in an outflow, it must overcome the local 
gravitational force of the galaxy's spheroid acting on the 
gas. The relative strengths of these two opposing forces can 
be estimated as follows. Consider a spherical optically thick 
cloud with radius R c , density n e , located a distance d from 
the nucleus. The gravitational force on the cloud is: 



7.5 Gas Motions in the EELR 



qrav 
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d 



(15) 



where m p is the proton mass, M(< d) is the spheroid mass 
within d and V c is the circular velocity corresponding to 
M(< d). 

The radiation force acting on this cloud is: 



F T ad — 



1 LphRc 

4c d 2 



(16) 



where L v h is the nuclear luminosity and assuming that all 
nuclear radiation illuminating the cloud is absorbed. The 
ram pressure force is: 



Fram = KPjV?R 2 



(17) 



where we used the definition of P r am from above. 

The ratio of gravitational force to radiation force on the 
cloud is: 



± gray 
Frad 
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(18) 



where we have used the cloud column density N e « 2n e R c . 
Significant acceleration by radiation pressure can only oc- 
cur if this ratio is < 1. Taking estimates of L p h and N e from 
above, and assuming a consta nt V c ~ 300 km s - 1 (based 
on the Ha rotation curve ffom lFerruit et aL | (|2004h as well 
as our own measurements of the southern NLR's [O III] ve- 
locities), this criterion implies a critical radius of d ~ 20 pc 
beyond which radiation pressure cannot effectively acceler- 
ate the gas. The inner plume deprojected size of I ~ 150 pc is 
much larger than this critical radius. This suggests that ra- 
diation pressure is not a significant acceleration mechanism 
for most of the resolved inner NLR. 

The ratio of gravitational force to ram pressure force is: 



Fran 



N e m p 



(19) 



This ratio is inversely proportional to nuclear distance. If 
Pram is roughly constant over the length of the jet, the ra- 
tio drops linearly with d, such that ram pressure begins to 
dominate over gravitational forces above a critical radius of 
about 20 pc. Unlike radiation pressure, ram pressure of the 
jet flow may be an important acceleration mechanism in the 
plume. 

A treatment of relativistic pressure in the same fash- 
ion is difficult, since it depends critically on physics of the 
boundary between the relativistic plasma and the ionised 
gas. Undoubtedly, gravity will have an effect on the displace- 
ment of any gaseous component in contact with relativistic 
plasma, but we forgo a discussion of relativistic pressure ef- 
fects, since the jet ram pressure is likely to dominate the 
acceleration of any gas that comes into contact with the jet 
flow. 



iFerruit et all |2004l ) argue that the asymmetry between the 
velocity profiles of the northern and southern EELR can be 
explained if the northern part is disturbed by the jet, while 
the south is in normal rotation in a thin disc. The low S/N 
of the line emission in our spectra on the scales of the EELR 
(nuclear radii > 1" ) does not allow a full exploration of this 
phenomenon. We can, however, confirm that an asymmetry 
does exist between the absolute peak velocity of the narrow 
emission in the NW and SE EELR along Slit A, measured at 
and beyond 1" (see Fig. 4). We may speculate on possible 
causes for the asymmetry. In our model for the NLR, the 
jet escapes the disk where most of the observed emission 
line gas is confined and therefore, direct interaction with 
the jet is unlikely to be the main cause of the asymmetry. 
This is supported by the fact that the jet is quite symmetric 
in its structure, which belies the basic symmetry of the gas 
density gradients that bend the jet. In this light, we propose 
two alternative hypotheses. 

A. ) We have shown that radiation pressure probably 
cannot drive a strong outflow beyond few tens of par- 
sees. However, it may still be able to influence the gravi- 
tational motion of AGN ionised gas on scales of the EELR. 
There is a strong difference in the mean ionisation between 
the northern and s outhern EELR (|Wilson fc BaldwirJli"985l : 
IFerruit et ail [2004). which implies a difference in the ionis- 
ing field strength and possibly the contribution of radiation 
pressure. A detailed study of the conditions of the emission 
line gas and any ionisation-kinematic trends in the EELR 
may test this idea, but this is beyond the scope of this work. 

B. ) The northern EELR has a pronounced dusty spi- 
ral arm that has no obvious counterpart in the south. Such 
structures have been proposed as paths along which gas 
streams into the nucleus from the outer parts of the g alaxy 
jMacieiewskill2004lFathi et al.ll2006l : lDavies et al.ll2009l ) due 
to spiral shocks. A comparison of the spiral pattern with the 
kinematic residuals left after accou nting for normal disc ro- 
tation , such as that presented by iGonzalez Delgado et al.l 
(2002), may provide evidence for this model. 



8 SUMMARY 

The circumnuclear environment of NGC 2110 consists of 
several components which each play a role in shaping the 
structure and dynamics of the region and the NLR. We ex- 
plore the outflow using well-resolved kinematics from STIS 
spectroscopy and constrain the influence of the circumnu- 
clear components using high resolution optical and radio 
imaging. Summarising our results: 

(i) A circumnuclear dusty disk plays a major part in shap- 
ing the appearance of the NLR. We propose a model that 
explains the angular offset between the axis of the inner 
radio jet and the EELR in terms of the biconical illumi- 
nation of this dusty gaseous disk. We also argue that the 
strong, symmetric bending of the jet is a consequence of the 
anisotropic mass and pressure distribution from the disk. 

(ii) A high surface brightness emission line structure, 
which we call the "inner plume" , is visible within the inner 
~ 100 pc of the nucleus. The kinematics of the inner plume 
reveal it to be accelerated, broadened and disturbed - the 
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site of a localised nuclear outflow. Its geometry indicates 
that it is not aligned along the radio jet, but possibly skirts 
the main jet flow. An ionisation analysis indicates that the 
inner plume is photo-ionised by the central active nucleus. 
Accounting for dust extinction in the disk, we conclude that 
the plume is possibly part of a extended symmetric outflow 
on either side of the nucleus. 

(iii) We explore the ionisation and kinematics of a bright 
compact nuclear knot, which has a size smaller than 20 pc 
and is barely resolved, even with HST. The low ionisation 
state of the knot and its high velocities, which are continuous 
with those of the inner plume, reveal that the nuclear outflow 
starts out on these small scales and then expands outwards. 
We speculate whether some of the kinematic structure in 
the knot may come from motion around the central black 
hole. 

(iv) Finally, we explore the dynamics of the outflow. Link- 
ing the jet bend to the pressure of the bulge atmosphere 
leads to constraints on the jet density and velocity. We sug- 
gest that this jet is significantly sub-relativistic and has suffi- 
cient ram pressure to drive the emission line outflow. Radia- 
tion pressure from the active nucleus or radio source pressure 
can both also account for the energetics of the inner plume, 
but the increasing influence of the gravitational potential of 
the galaxy with nuclear distance restricts the importance of 
radiation pressure in accelerating gas in the resolved inner 
plume. 
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